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OBJECT 


To  determine 
exploding  body  on 


the  influence  ef  the  surfaoe  oontour  of  an 
fragment  distribution. 


ABSTRACT 


The  effeot  of  the  surfaoe  contour  of  an  aerial  exploding  body 
upon  fragment  distribution  is  derived  for  a surfaoe  of  revolution 
of  arbitrary  shape* 


CONCLUSION 


The  density  of  the  fragment  distribution  at  any  position  on  a 
target  area  can  be  determined  for  any  given  symmetrical  surfaoe 
contour  of  a missile.  Generally,  this  distribution  -will  be  a vari- 
able. 


The  primary  interest  in  desiring  an  analytical  expression  for 
fragment  distribution  is  to  determine  its  effect  on  lethal  area. 
Naturally,  the  most  desirable  can torn*  would  be  that  which  gives  the 
may! puna  lethal  area  for  a given  number  of  fragments. 

From  the  equations  derived  in  this  report,  it  is  possible  to 
determine  the  fragment  distribution  of  a missile,  given  the  equation 
or  coordinates  of  the  missile  surfaoe.  If  the  fragment  distribution 
is  known,  the  lethal  area  for  the  particular  taotioal  application 
intended  for  the  missile  oan  be  calculated  by  using  the  information 
furnished  in  various  Ballistio  Research  Laboratories  (BRL)  reports. 
In  many  oases,  the  computations  owild  be  set  up  in  tabular  form,  so 
that  a computer  oould  carry  out  the  calculations  without  understand- 
ing the  physios  of  the  problem. 
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SYMBOLS 


Variable  distance  along  x axis 
Variable  distance  along  y axis 
Distance  from  origin  to  element  ds 
Angle  between  r and  the  y axis 

Distance  from  x axis  to  the  center  of  the  radius  of  curvature 
for  any  element  ds 

Radius  of  curvature  for  rny  element  ds 
Angle  between  p and  the  y axis 

Distance  from  x axis  to  the  lowest  point  of  the  body 
Vertical  distance  from  the  lowest  point  of  the  body  to  the 
target  surface  T 
A / h 

An  inf initesimal  arc  length 

Distance  from  the  y axis  to  the  x- intercept  of  the  curve  s 
Distance  from  the  y axis  to  any  arbitrary  point  an  the 
target  surface 

(-ay1  / yy1  / x) 

Diameter  of  the  desired  target  area 

Distance  from  the  center  of  curvature  of  any  element  ds  to 
point  a on  the  target  surfaoe  (see  Fig  1) 

Surface  area  of  the  exploding  body 
Projected  area  of  the  fragment 

Number  of  fragments  per  unit  of  the  exploding  body 
Density  of  the  fragment  distribution  - number  of  fragments 
per  unit  area  impinging  on  the  target  surfaoe 
Weight  of  a fragment 


dy/dx 

d^y/dx^ 

(1  / (y1)2  “ ay11  / yy11) 


R - Distances  fragment  travels  from  missile  surfaoe  to  target 
T - garget  surface 
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INTRODUCTION 

1.  The  expeoted  number  of  fragment  hits  on  a target  depends 
on  the  fragment  density  at  the  target  and  the  projected  area  °ff®^e 
by  the  target.  Whether  the  target  is  incapacitated  or  not  dep< ends 
on  the  mass,  velocity,  and  shape  of  the  fragments,  in  addition  to 
the  number  and  location  of  hits.  This  report  will  be  concerned  with 
only  one  of  these  factors,  the  fragment  distritaticm.  Since  the 
most  important  factor  in  controlling  fragment  distribution  is  t e 
surface  contour  of  the  exploding  body,  the  analytical  approach  to 
the  effect  of  the  surface  contour  is  considered,  restricted  by  tne 
following  assumptions: 


a.  A constant  number  of  fragments  per  unit  surface 

area  of  the  exploding  body. 

b.  The  fragments  leave  normal  to  ihe  surface  of 

the  body 

c.  The  velocities  of  the  fragments  are  of  such  a 

high  value  that  the  influence  of  gravity 
can  be  neglected 

d.  The  exploding  body  is  considered  stationary; 

that  is,  the  body  velooity  is  considered 
to  be  zero. 


In  order  to  obtain  a simple  mathematical  solution  to  the  problem, 
certain  assumptions  must  be  made.  It  is  believed  that  the  above 
assumptions  arc  justifiable. 

2.  The  first  assumption,  of  course,  can  be  controlled  in  the 
fabrication  of  the  projectile. 

3.  The  second  ascumption  is  oertainly  the  most  debatable. 

The  location  of  the  point  of  detonation  will  xmdoubtedly  affect  the 
normal  radiation  to  some  extent.  If  this  influence  is  constant  for 
a given  design,  however,  as  it  should  be,  the  fragment  distribution, 
considering  normal  departure,  should  be  susceptible  to  mathematical 
correction.  All  other  influences  on  normal  departure  from  a static 
body  should  be  of  a random  nature  and  hence  should  fluctuate  about 
the  normal  as  an  average. 

I;.  The  third  assumption  restricts  the  fragment  paths  to 
straight  lines. 


3 

\ 

CONFIDENTIAL 


CONFIDENTIAL 


5.  The  fourth  assumption  is  related  to  the  second  assumption 
in  that  any  velocity  of  the  tody  will  affect  the  normal  radiation, 
of  the  fragments*  In  many  cases,  however,  the  velocity  of  the 
body  will  have  little  effect  and  can  be  neglected. 

6.  The  mthematics  of  this  paper  were  completed  before  the 
author  realized  that  parts  of  it  paralleled  analytical  work  aocomr 
plished  at  the  Ballistic  Research  Laboratories.  The  surface  dis- 
cussed in  paragraphs  13-16  of  this  report  is  known  as  the  Kent- 
TTitchcock  contour.* 

7.  m addition,  H.  K.  Weiss  **  of  the  Ballistic  Research  Labo- 
ratories has  derived  the  optimum  contour  for  certain  tactical  sit- 
uations* 


* Kent,  R.  H. , "The  Shape  of  a fragmentation  Bomb 

to  Produce  Uniform  Fragment  Densities 
on  trie  Ground,"  BRL  Report  No.  7&2 

**  Weiss,  Herbert  K.,  "Optimum  Angular  Fragment  Dis- 
tributions for  Air-Ground  Warheads," 
BRL  Report  829 
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DISCUSSION  OF  RESULTS 


General  Solution  for  Axis  Perpendicular  to  Target  Surface 
8.  General  solution: 


Consider  that  a missile  of  arbitrary  shape,  symmetrical  about  its 
axis,  is  detonated  with  its  axis  perpendicular  to  a plane  target 
surface  (see  Fig  1).  The  density  of  the  fragment  distribution  on 
the  target  surface  T for  any  position  ' 3 given  by 


K -/ldAg 


dArj 

(i) 

where 

dAa  m 2 TT'xds 

(2) 

and 

dAq  r 2<hn 

(3) 

Substituting  Equations  2 and  J in  Equation  1 gives 

TT  - /*!  *^8 

(b) 
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..  . J.  n and  dn  &s  functions  of  x,  y, 
9,  It  now  remains  no  find  ds,  r\  , H 

yl  and  y11 

From  Figure  1.  it  can  he  seen  that 

n s (»/p)  rtn  *-((»■*»«-'  sin0° 

and  that  (R  /f>)  - cos  Tg 

Suhstituting  Equation  6 in  Equation  ^ there  results 
q - (a  - h)  tan  6 - p sinV  r sinOC 

Ic  is  also  noted  that,  * * * oos0t  “ l°°°S  ® 

Suhstituting  Equation  8 in  Equation  7.  there  results 
n.atanfl-r  cosflC  tan  6 / r since 


(5) 

(6) 

(7) 

(6) 

(9) 


From  elementary  calculus  the  expression  for  the  radius  of  curvature 
in  rectangular  coordinates  is  given  y 


P = - (i 


(10) 


Also  from  Figure  1,  it  is  that 


-dy  “ tan  © 

3x 


Therefore# 


3in  9 ’ ~(T  7 &rrjrfir 


(11) 

(12) 


10.  By  using  Equations  10,  11,  and  12,  and  since  y = * coso<L 

• wmiation  9 can  he  written  as 

and  x ■ r sinOC  , Equation  7 


- ay1  / yy1  / x 


(13) 
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The  following  shows  how  df}  varies  with  x,  y and  yl 

to  /<^y  / dy1 

d y 

dr}  = £n  d*  / aoy1^  / dr^yixdx 


(14) 


(15) 


r 


(16) 


(17) 


Considering  Equation  13,  Equation  15  becomes 

dq  c (1  / (y1)2  “ ay 11  / yyll)  dx 

Furthermore,  ds  can  be  expressed  as 

ds  m (1  / (y1)2)1/2  dx 

Finally,  substituting  Equations  13*  16,  and  17  in  Equation  i+ 

there  results  K mS\  * (1  / fy*  )2^  ;2  __  (18) 

(x  - (a-y)  ylj  Q-/  (yl)d  - (a-y)yilj 

From  Equations  13  and  18,  the  density  of  the  fragment  distribution 
on  target  T can  be  obtained  at  any  position  ty  for  any  given  sur- 
face contour  of  the  exploding  body. 

11.  The  expression  for  the  fragment  distribution  (Kn)  normal 
to  its  path  R at  point  p can  be  expressed  simply  as: 


K 


cos  © 


or  since  cos  © - 


rr  / (yDSJl/6 


and  considering  Equation  18 
whereas  ■ -ay^  / yy1  / x 


Kn  s six  (1  / (y1)2) 

V *)■  XT  777)?  -ayll  / "yyll) 


(19) 

(20) 

(21) 

(22) 
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Simplification  of  general  Soluti on 

12.  Although  the  solution  obtained  in  the  previous  aeotion 

forrsr?ican?  °rrt-iD 

by  atS£g  P aBd  nearly  " ,00urate-  Thi»  *». 

R 

r fox\ 


therefore  q s (KJ  sin  « • a tan  « « ayf 
and  a 

ar?  ® “ay11dx 

As  a result.  Equation  i+  can  be  written 
X = y«yx(l  / (yl)2)l/g 


also  Kn  B / (y^)2) 

a Jrly  (27) 

Where  in  both  cases  q „ >ayl  (2Q) 

This  solution  mould  no  longer  be  * 

1*iohf’ oreafc«r  SUrfa°e 

Jgaot  Solatia.,  far  tha  Sarfaaa  that  Produce.  Canat«t  Frag.aat  Dan.ity 

S3.^  pro-~ 

the  target  area  T*.  The  ideal  approach  +Q  +M  « radiating  upon 
formally  solve  the  nonlinear  differential  E wati STls' bevto 
a constant.  This  equation  of  course  is  udm^^y^L.*******6  K as 
constant.  There  should,  home^T  -oSS1o\1t^\\r^°ld  7 * 

* See  Appendix  A 
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result  in  a curved  surface,  the  surface  in  *ioh  iw  are  interested. 
This  can  be  obtained  by  solving  the  approximate  Equation  2~,  tr.a 
ingK  as  a constant,  and  keeping  in  mind  the  restrictions  imposed 
upL  this  equation*4.  The  order  of  approximation  of  the  solution 
of  Equation  26  can  be  obtained  by  substi tutang  i*®  d®riv®f  surface 
in  *CIZ  18.  *us  obtaining  the  exact  density  of  the  ^gmeut 
distribution  produced  by  this  surfaoe.  Equation  26  can  be  written 

K a2  p dp  ■ x (1  / p2)1/2 
"dx 


where  p 


■4 


Transforming  ka2  f(l  / P2)  1^2  (2 

Spif  J 

Upon  integrating  ka2  (1  / P2)1^2 


(2pdp) 


x2  / c 
2 


ll*.  Wow  consider  the  boundary  condition 


- 0 . Then ' 


Substituting  for  C in  Equation  31  gives 


C • ka^ 

p*r 


ka2  (l/p2)l/2  --  # 


2 

Squaring  the  above  equation  on  both  sides  and  solving  for  p, 

there  results  o , 0.0  , 

p2e/r2  (x2  /2ka2)2  -1  , 

yr~  (3U) 

Uk?  ah 

Simplifying  and  solving  for  p gives 

p - - (x2  / 2)l/2  x (35) 


It  is  of  interest  to  note  that  the  solution  y a constant  no 
longer  satisfies  this  equation. 
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^ero  the  negative  sign  ia  chosen  outside  the  radioal  to  conform 
with  -the  chosen  coordinate  system  and 


Ska2  e A. 


remembering  that  p c dy/dx,  Equation  35  becomes 
V a - f (x2  / 2)1/2  x dx 

J 7 


Upon  integrating 


A, 

y z -A,  (X2  / 2)3/2  / C1 

3 ? 


(36) 

(37) 


or 


15-  Consider,  now,  the  boundary  condition 

y)  r 0 where  A.  is  now  considered©, 
x -»  0 

= "A,  8 / Cj 

3 

s : 0 / A,  8 

7 


(38) 

(39) 


!‘r°m  E^uat,ian  to  into  Equation  39  leads  to  the 
deseed  surface  (Equation  1*0).  which  will  produce  a extent  den sitv 
of  fragmmta,  subjected  to  the  assumptions  made  in  pa^^h  £? * 


y r -A, 


((*£/2)3/s  -23/2)  (See  Fig  3) 
3 A2 


(J|0) 


where 

and 


A, . ^ 


q : a (x2  / 2) i/s  x 

1 __  _ 


Af 


(to) 


n 160  If/or  Pra°tical  construction  reasons  (for  example  to 
allow  space  for  a fuse  well),  a diameter  (d)  at  the  ap^S  the 
surface  must  be  inoperative,  it  would  be  desirable  to^use  the 

oX^  , lAm  \ =D/2  ^en  solving  for  the  constant  C,  in 

order  to  cut  down  on  the  dead  space  on  the  target  surface.  1 
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Using  this  boundary  condition  in  conjunction  with  Equation  Jl,  there 

results  9 

C - 4?  - 

2 * “ff  (i£) 

Consequently  A^  (1  / p*1)'1'/2  s x2  / A?  - D2  (1+3 ) 

t 5"  r*  T 


(*2  / l ~D2a2-1 

A*  43* 


Equation  14*  is  the  solution  in  integral  form  for  this  particular 
case.  Given  constants  D and  A,  "tile  desired  surface  is  easily 
obtained  by  numerical  integration. 


Fragment  Distribution  for  a Spherical  Missile 

17.  The  equation  for  a oircle  in  polar  coordinates  is 
y = cos  c 

x b sin  © (1*5) 

where  r is  a constant.  When  revolved  about  a diameter  a spherical 
surface  is  generated. 


y1  = 

ay- 

d©  s 

sin  © ■ tan  0 

"aar 

cos  © 

(1+6) 

■H  - 

38 

dy1 

33“  - 

«>sec2©  - -sec3© 

Off) 

“cET 

38 

r cos  © r 

Substituting  the  valixes  for  x,  y^-,  and  yl  I in  Equation  26,  there 
results 

K r y»fr2  cos?  © ( 1 o\ 


: SI  r* 


CONFIDENT1AL 


•*^*s 


it  • 


siteiitS® 


CONFIDENTIAL 


Letting 


-*4  2 A 


then  K s /\cos^6 

and  expressed  as  a dimensionless  ratio 


K a co  s3©  (See  Fig  h) 

A 

Ceneral  Solution  for  Ax is  Inclined.  to  Target  Surface 


18.  The  solution  obtained  in  paragraphs  8 to  11  is  restricted 
in  that  the  axis  of  the  missile  rust  ho  perpendicular  to  the  target 
surface.  Obviously,  this  is  not  always  the  case,  and  it  becomes  of 
interest  to  find  out  how  the  fragment  distribution  varies  v/hem  the 
missile  axis  is  inclined  to  the  target  surface. 
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19.  Consider  Figure  2.  The  plane  (eo'd)  is  perpendicular  to 
line  (oo’)»  It  is  desired  to  find  the  fragment  density  distribution 
cn  plane  (go'd)  ■which  is  inclined  to  line  (oo1  ).  The  density  K on 
plane  (eo'd)  at  any  point  b is  given  by  Equation  18,  or  by  Equation 
26  if  the  simplified  solution  is  justified.  The  fragment  density 
normal  to  R at  point  a is  therefore 


IL  = K_ 

na  cos< 


However,  (ob)  can  be  expressed  as 


and  (oa)  as 


(oa)  z a sin 


coax 


sin  ( TT  - e -OC)  cos  (©  /<X)  (53) 

“2 


Kn  can  therefore  be  expressed  as 
a. 


Kn  = K cos2  (e  7/oC) 


cos^doos 


and  the  fragment  density  on  nlane  (go'd)  at  point  a as 

K_  " K>oos2(«/pC)  (55) 

3 5 

oos^cos  OC 

It  is  more  convenient  to  express  the  angle  Ot  as  a function  of  , 
the  angle  the  missile  axis  (00* ) makes  with  the  normal  to  the 
target  surfaoe  (go'd),  and  the  angle  between  planes  (ogo1 ) and 
(oao1).  This  can  be  done  as  follows 

sin  s bh 

To  (56) 

But  be  • a tan  © cos  (57) 

and  bh  - a tan  © sin^  (58) 
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Consequently 


sin^C  - sin  0 cos  fi 


tribution  S»^dfr  J°  loca^e  the  position  of  our  new  fragment  dis- 

{I  ? v tennS  0f  (o'h)  mst  also  be  determined.  From 
figure  2,  it  is  observed  thati 

— (o'b)  f6o) 


(o*a)  - 

cos  © 


(o*b) 
cos  (P 


7 ?f  hbf  esneral  solution  in  paragraphs  8 to  11,  so 
by  letting  (o»a)  ■ * 


there  results 


n'  r noosd 

cos  (©  /ot ) 


811(1  locatian  oan  b«  determined  on  the  inclined 
Jjf16  ?y  ^ing  B(3uatiQas  55,  59,  and  62,  as  long  as  the 
distribution  is  known  on  the  plane  normal  to  the  missile  axis. 

the  WbbS^°^Jl!e+Stri^g  of  a fragment  is  a function  of 

also  of  £b«r^  graveled,  the  length  of  the  new  fragment  path  is 
also  of  interest.  It  oan  be  expressed  as 

(°a)  s a ocb PC  - h cosctC 

cos  /pc  ) cos^cos(e  /oC) 

bo^b611^  °f  ?e  !lag?ent  di^ded  by  the  length  of  the  path 

to  the  normal  target  plane  is  therefore  * 

(°a)  z Q08  3 oosflC.  (Ah) 

ToBJ  cos  (©  / *C) 

ofP  f2*.  Rations  55,  59,  and  62  must  be  used  to  determine  the 
£2? ' “iS  °f  the  3haM  d*ri'r*d  10  17. 

Ka  - g COs^  (*  /gC ) 
cos^  ^ OOo^'^C 

sin^,;  sin  ^ oos^J 

n!  * q cos  e 


where  for  this  particular  problem  K is  a constant. 
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23#  P or  a given  angle  $ j end  ooula  be  plotted  as 

functions  of  © for  various  angles  A . If  desired,  a set  of  tables 
or  graphs  could  be  drawn  up  in  tills  manner.  It  would  then  be  a rather 
simple  process  to  compute  the  lethal  area  as  a function  of  the  angle 
of  inclination. 


Liolosuresi 

1-2  Figures  3 and  1* 

5-7  Appendioes  A,  B,  and  C 
(Appendix  B contains 
Tables  1 and  2) 
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*'J*£p*  • Acos*6  • 

zv  * frsf/umt  ares  on  m/ss/M 
surface 

r * radius  of  missile 
S • ho/fhf  of  burst 
X*  fray/ unit  ares  on  tarfml 
surf aco 


fragment  distribution  of  a spherical  shell 

, UPON  A PLANE  TARGET  SURFACE 
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APPENDIX  A 

SimF ACE  THAI  PRODUCES  OPTIMUM  LETHAL  AREA  UNDER  CERTAIN  ASSUMPTIONS 
(This  problem  was  solved  by  Mr.  H.  K.  Reiss'*) 


The  assumptions  stated  in  the  "Introduction"  are  assumed  to  hold, 
and  in  addition  the  axis  of  the  missile  is  considered  perpendicular 
to  the  target  surface.  The  expected  number  of  hits  on  a target  can 
be  expressed  as: 

Eh  = a© 

r'd  (la) 

where  o fragments  per  unit  area  per  unit  distance 

from  the  burst  point,  at  angle  © 

r 2 distance  the  fragments  travel  to  ihe  target 

A©  e exposed  area  of  idle  target 


The  expected  number  of  disabling  hits  is 

Ek  e Eh  Phk  (2&) 

where  Phk  is  the  conditional  probability  that  a hit  will  disable  the 
target.  This  is  usually  a function  of  the  mass,  velocity,  and  the 
presented  area  of  a fragment.  For  mathematical  simplicity,  P^y  will 

bo  treated  as  a constant.  This  means  that  all  ihe  fragments  are 
considered  to  have  ihe  same  disabling  potential. 


By  using  the  Poisson  approximation,  the  probability  of  disabling 


ihe  target  is 


-Bj 

Pk  = 1-e 

#o  *o 


(3«0 


The  number  of  disabled  targets  is  ^ ^ L*rpk  dxdy3  ""here  O' is  the 


number  of  targets  per  unit  area.  If  O"  is  considered  oonstant,  it 
can  be  brought  outside  the  integral.  The  integral  term  is  then 
defined  as  the  lethal  area  (A^). 


09  tO 


al  « j"  J"  pk  tody 


Weiss,  H.  K.e  op0  cit. 


-•©  -to 
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If  the  lethal  area  is  multiplied  by  the  density  of  the  targets,  there- 
fore, the  number  of  targets  disabled  will  be  obtained.  Equation  U&  can 
be  expressed  in  polar  coordinates  as 

m 2f>a^^J"  (l-e~*k)  sin  0 seo3  © d 0 (5a) 

a a height  of  burst 

0 - angle  between  missile  axis  and  fragment  path 


wnere 


It  is  convenient  at  this  point  to  make  the  following  change  of 
variables 

IT  — cos  0 

dW  a -sin  0 d 0 


and 

Equation  5a  becomes 

Furthermore 

However 

Therefore 


Al  = -2na2  f (!-•  ”Ek)  r5  dw  (6a) 
= Eh  phk  = M(9>  A -iiL  pMe  (7a) 

r r a/oos  0 b aw-1 
Ek  = M(e)  A (0)  H2  Phk  (8a) 


If  it  is  also  oansldered  that  the  exposed  area  varies  as  the  00s  0, 
as  would  a flat  target,  Iq  can  be  expressed  as 


Hence 


and 


Ae  : A 00s  0 s Aw 
E]r  E \ Py 


"3*" 


Al  = -2ira 


J (1- 


hk 


rhk 


) w dw 


(9a) 

(10a) 

(11a) 


A sail  1 increment  of  area  at  some  position  (W),  over  a band  width 

(6  *)  is  , 

(c5Al)  - -2  IT  a2  (1-e  hk)Tf5  <5  W (12a) 
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If  M is  now  increased  at  the  point  (W),  over  the  band  width  (<5W) 
the  increase  in  (t/AL).  is 

r phkl  . 

6 ( 6 AJ,  e I 2 *rr  A e — — dw  <5m  ( ’ 3a) 

^ a^ 

but  since  the  total  number  of  fragments  must  not  increase, 

M must  be  subtracted  at  some  other  point  (Wg)*  At  that  point 

j -Hr  AW^  pJ^WcfM  (Ilia) 


d(cTAL) 


2*A  phk  e 


The  necessary  condition  for  to  be  the  surface  that  produces 
maximum  lethal  area  is 

6(6  al)1  -6(6  al)2  = o (15a) 


Therefore 


21CA  . 


-IVAlfJ  Pht 


s C ■ constant  (16a) 


P-.  r Ln  r 2 'Tt* A P. 


(17a) 


s f a2  Ln  ( 2 Tf  A.  P, , )1 

K7  -TT^I 


a2  Ln(2irAP,.)  s A* 


\ e A»W"3 


(18a) 


(19a) 


- A*  sec^  © 


(20a) 
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Consider  now  the  Eent-Hitchcock  contour 
requirements  "tJiat  K ~ constant  where 


K 


■ ytfxds 


or 


But 


K - Ife  sin  9 d 9 
rj  s a tan  0 


It  fulfilled  'the 

(21a) 

(22a 


drp 


a soc^  8 d 9 


Therefore 


K r sin  9 d 8 

„2  tan  0 sec2  8 d 8 


Ka^  sec3  0 


(2?a) 


(2l4a) 


or  can  he  expressed  as 

Mg  . 

Equation  2lja  has  the  form  of  Equation  20a.  Therefore,  the  Kent- 
Eitchcock  contour  is  the  contour  that  produces  optinum  lethal  area 
under  the  restrictions  of  this  particular  problem. 
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APPENDIX  B 


APPLICATION  OF  ANALYSIS  TO  AN  AERIAL  EXPLODING  BODY 


1 . 


i 

j 

t 

it 


Problem 

find  the  proper  surfaoe  oontour  of  an  aerial  exploding 
body  that  -will  produce  a fragment  distribution  of  constant  density 
over  the  target  surface  T under  the  following  restrictions: 

1.  Height  (h)  of  the  exploding  body  above  target  - 20  ft 

2.  Desired  density  of  the  fragment  distribution  (k)  over  the 

target  T - 1-67  fragments/ffc2 


3.  Desired  diameter  (D)  of  target  r h®  ft 
I;.  Fragment  rreight  z 2.5  grains 


5. 

O. 


p 

Surface  area  of  fragment  kf  e -01175  in 
Fragment  distribution  (/t)  on  body  surfaoe  r 
1 r 85-106  fragment s/in2 

.01175' 


The  surface  which  satisfies  tlrs  problem  is  that  of  Equation  1*0. 
For  simplicity,  however,  let  a = h and  therefore  A e 
Equation  I4.O  would  then  reduce  to 

(lb) 


y = -A  [(x2  / 2)3/2  - (2)^2] 


In  order  to  establish  the  boundary  of  the  body  in  the  x direction, 
an  additional  equation  concerning  the  surface  area  of  the  body  is 
needed.  This  is  expressed  math era. tically  as 


As 


r 2^*  tt(l  / (y1)2)1/2  xdx 


where  from  Equation  33  it  can  be  shown  iiiat 

(1  / (y1)2)1/2  s (x2  / 1) 

IT 
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(3b) 
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Equation  2b  can  therefor©  be  expressed  as 
As  a 2 


I ° % / x) 

c 


Upon  integrating 


As  - 


(xcA  / XO2) 


2A 


(Lb) 

(5b) 


Equation  5b  oan  be  expressed  as 


XoL  / 2A2  xo2  - 2A2As  r 0 (6b) 

•It* 

■where  A2  z 2ka2  b 2(1.67)(20)2  c 15*698  in2 

7T-  — 


and  in  this  particular  problem 


As 


(1.67)(. 01175)  e 2L.658  in2 


Substituting  these  values  in  Equation  6b  and  solving  for  xQ 
by  using  the  quadratic  formula,  it  is  found  that 


Xo  - 2*550  in. 

By  substitution  of  the  value  of  A in  Equation  1+0,  the  desired  surface 
becomes 

y m -I.52O0  [ (.06570x2  / 2)5//2  -2.8281+]  (7b) 

from  uhicii  y^  r -.25237  x (•06J10j£  / 2 )^/2  (8b) 

and  y^  — -.01608  x2  (.06370X2  / 2)“^/2  -.25237  ( .O637OX2 

/ 2)1/2  (9b) 


In  order  to  find  the  exact  daisity  of  fragment  distribution 
on  the  target  surface  produo ed  by  Equation  7b,  it  is  only  necessary 
to  substitute  the  above  expressions  in  Equation  18.  The  results  of 
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this  substitution  are  presented  in  Table  1*  It  can  be  seen  that  'he 
density  of  the  fragment  distribution  is  approximately  2#  to  1$  below 
the  desired  density  of  the  fragment  distribution.  It  is  also  seen 
that  the  target  area  is  3/£  larger  than  specified.  This  result  sug- 
gests that  the  derived,  surface  equation  (lb)  can  be  corrected 
to  give  improved  results*  It  is  proposed  that  the  density  of  the 
fragment  distribution  be  increased  3%  by  specifying  that  K s 1.03K 

Then  z 2 (1*0?  ka2  m (1.03)(15*698)  » 16. 169 

A r lu  021 

Hence,  the  corrected  equation  for  y becomes 

y *-1.34030  (.06185X2  / 2)3/2  - 2. 8281; 
and  yl  - -.2ijB69x  (.06185*^  / 2)V^ 

y11  r “*01538x2  (.06185x2  / 2)“1/2  -.24869 

(.06185x2  / 2 )x/2  (12b) 

Table  2 gives  the  improved  results.  The  density  of  the  frag- 
ment distribution  now  varies  from  / l/2/S  to  - 1%,  while  the  target 
area  is  in  error  by  less  than  l/2?& 


(10b) 

(Hb) 
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APPENDIX  C 


pro liE  2rfot^ft™risEittrofc  COntOT 

area Ts  t^  V ^r  ho!r  muok  6«M«r  «*,  lethal 
tarjet  radius  (L) 


The  distribution  of  fragments  from  a spher 


e 11 


X 

s 


K 


o 


COS- 


e 


(lc) 


The  number  of  fragments  over  the  target  surface  or  radius  L is 
^ "L 


r 2*rtE0  1 cos3  €>  rdr  m 


f 


r? 


rdr 


(r2  / L2)3/2 


(2c) 


the  iMSLS  *“=*  " "*«  L Proa 


Kkh  = *"lch  *n*L£ 


(3c) 
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If  if 


-v<  K ^ 


•^a»  % 


^ «*  O, 

-v.  as 
• > Sr  c 

” c ; «-  xs  m 


■ Fragment  Distribution  Produced 
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To  make  a comparison  of  lethal  area,  is  assumed  to  be  equal 
to  N . By  equating  these  two  expressions,  it  is  found  that 


K0 

S Z 1*706Ky1 

(2-  y(2) 

(ho) 

Lethal 

area  can. 

be  expressed  as  (Equation  6a) 
- -2  <l-a2  f (l-e"Ek)ir3dw 

(5c) 

where 

Be  C 

(6c) 

Consider  the  Aj  ; A cos  6 = AFT.  The  lethal  area  for  a 

sphere  is 

AL/s 

- 2fTA2 

J (l-e-1*?06  A^khPhk^  *~5dTr 

(7o  ) 

since 

K > 

K0w5  - 1.706E1&W3 

where 

A^kh 

P^  r constant 

The  lethal  area  for  tiie  Kent-Hitohcock  contour  can  be  expressed  as 
Al  | - 2 era2  J"  ( 1 - e-^kh^k^)  ir3dw  (8c) 


Assume  AK^P^  ; 2,5 

(This  is  a fair  value  for  a 

sample  problem) 

Then  A^  i f 1 

[ (i*-88*)  ir-3*. 

(9o) 

2ora2  ' 8 

20a2  'ith 

(1-e  -2.5F)  w~3dw 

(10c) 
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Further  assuming  that  the  limit  of  the  target  area  corresponds 
to  © 6 k5°t  then  the  limits  on  the  above  integration  are  1 to  «707 
since  W - cos  ©.  Carrying  out  the  above  integration  numerically, 
there  results 

AL/kh  ® l»l66  (He) 

AL/s 


That  is,  the  Kent-Hitchcock  contour  gives  a lethal  area  that  is 
l6»6%  greater  than  a sphere,  assuming  that  both  surfaces  contain  the 
same  number  of  fragments. 
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